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An Accurate Theoretical Model for the Thin-Wire 
Circular Half-Loop Antenna 


Guangping Zhou, Student Member, IEEE, and Glenn S. Smith, Fellow, IEEE 


Abstract—The conventional Fourier series analysis for the 
thin-wire circular transmitting loop, or its image equivalent the 
half-loop, uses a delta-function generator for excitation. This 
method of excitation introduces two problems: it does not 
correspond to any realizable method of feeding the antenna so 
an accurate comparison with measurement is not possible, and it 
produces a divergent series for the input susceptance. To over- 
come these problems, a new theoretical model is used for the 
antenna: a half-loop driven through an image plane by a coaxial 
transmission line, with a transverse electromagnetic (TEM) mode 
assumed in the aperture of the coaxial line. This model is solved 
in a manner that preserves the simplicity of the original Fourier 
series analysis; all coefficients are obtained as closed form ex- 
pressions. Input admittances calculated from this new theoreti- 
cal model are in excellent agreement with accurate measure- 
ments. 


I. INTRODUCTION 


HE circular-loop antenna shown in Fig. 1 has radius b 

and is constructed from a perfectly conducting wire of 
radius a;. When the wire is thin, a;<b, 8,a; <1, the 
time-harmonic (e/“’) surface current on the loop is predomi- 
nantly in the azimuthal direction, ¢, and is conveniently 
expressed as a Fourier series: 


+ I(¢). 
Js > 27a; ¢ 
I($) =In +2 >. I, cos(ng). (1) 
n=1 


Apparently the Fourier series was first used for the thin- 
wire circular loop by Pocklington in 1897 to study the 
scattering of an incident plane wave by a closed loop [1]. 
Later Hallén [2], Storer [3], and Wu [4] applied the Fourier 
series to the thin-wire circular transmitting loop; the excita- 
tion being a time-harmonic incident electric field concen- 
trated at @ = 0 on the surface of the loop: 


E4(¢) = V,8(¢)/6, (2) 


the so-called ‘‘delta-function generator.’’ A great accom- 
plishment of this theory is that the Fourier series coefficients, 
I,, for the transmitting loop are available as closed form 
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Fig. 1. Circular transmitting loop excited by delta-function generator. 


expressions [4]-[7]. In the formulation of Wu, they are 


pec 0,1,2 3 
nn 3 a, ’ n= sto ms ( ) 
with 
6,0 n? 
a, = 2 (Kua + K,_1) = ,b *” (4) 
and 
1 a; a; 3 
Ky = <5 Vin? = (800) | x6] Vo? = (8,0) | 
1 1/2 
oer 2 — R2p2 
+— fy + in [4(n? - 6302)!" } 
2 nu} 1 
me fry Um+1 
1 26,6 
-> [9,,(x) + JJ2,(x)] ax. (5) 
2 x=0 


Here %, and %) are the modified Bessel functions of the 
first and second kinds and order zero, 0,, is the 
Lommel-Weber function of order 2”, J,, is the Bessel 
function of the first kind and order 27, y = 0.57722 -+> is 
Euler’s constant, and 6, = 27/A, = w/c, &, = Vbo/Eo- 

The Fourier series analysis for the thin-wire circular loop 
antenna has been extended to treat several other antenna 
structures, these include: the loop with discrete or continuous 
impedance loading [8], [9], the loop with a spherical ferrite 
core [10], loops with concentric and eccentric spherical 
insulation in a dissipative medium [11]-[13], and the loop 
over a material half-space [14]-[17]. In all cases the delta- 
function generator is the excitation for the transmitting loop. 
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Fig. 2. Comparisons of input admittances as function of the number of 


terms in Fourier series for two antennas with different generators. 8,b = 
1.0, 2 = 10.0, a, /a; = 2.30 (Zp = SO Q). 


10000 


The use of the delta-function generator (2) leads to a 
well-known problem, the imaginary part of the series for the 
current (1) is divergent at the drive point: Im[/(0)] > 0 as 
n— oo, In terms of the input admittance 


Y=G+jB=1(0)/V, = [+23 Ae (6) 
n=1 


this means that the series for the input susceptance B is 
divergent. This divergence is clearly shown in Fig. 2, where 
the admittance (G, B) is graphed as a function of the number 
of terms in the series for a loop with 6,b = 1.0 and Q= 
2In(27b/a;) = 10.0. For practical calculations the series 
for the admittance (6) is usually truncated at 20 terms, and 
reasonable agreement is obtained with experimental values. 
Clearly this problem with the convergence detracts from the 
otherwise elegant theory for the thin-wire circular loop an- 
tenna. 

To eliminate the divergent susceptance, Inagaki ef al. 
replaced the delta-function generator by a pulse-function 
generator: a combination of two Heaviside unit-step functions 
that maintains a uniform electric field across a small gap 
in the conductor [18]. The problem with this model is that, 
like the delta-function generator, it does not correspond to 
any realizable method of excitation for the loop. One’s initial 
thought may be that this model corresponds to a circular loop 
driven by a parallel-wire transmission line. However, the 
electric field between the conductors of a parallel-wire line is 
not uniform; and at the end of the line, the electric field is not 
rotationally symmetric about the axis of the loop conductor. 

The half-loop antenna consists of a half circle of wire 
mounted vertically over a conducting image plane. The half 
loop is the image equivalent of the full loop. It can be 
analyzed using the procedure presented above when the 
exciting electric field is a delta-function or pulse-function 
generator applied between one end of the loop and the image 
plane. The current and input admittance of the half-loop 
antenna are then precisely twice those for the full-loop an- 
tenna (1), (6). 

In this paper the aforementioned problems are eliminated 
by using a more realistic model for the excitation of the 
half-loop antenna: one that corresponds more closely to a 
practical method of excitation and does not lead to a diver- 
gent series for the susceptance. The Fourier series coeffi- 
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(a) Half-loop driven through image plane from coaxial transmission 
line. (b) Equivalent full loop excited by magnetic frill. 


Fig. 3. 


cients for the current are given by closed-form expressions 
that are no more complicated than (3)-(5). Input admittances 
calculated from the new theoretical model are in excellent 
agreement with accurate measurements. 


Il. MopEL FoR THE ANTENNA 


The new modei for the antenna is a half-loop driven 
through a perfectly conducting image plane from a coaxial 
transmission line, as shown in Fig. 3(a). The electric field in 
the coaxial aperture is assumed to be that of a transverse 
electromagnetic (TEM) mode: 


E,(e) = V,/[ein(a,/a;)], (7) 


where p is the radial coordinate on the image plane, and a; 
and a, are the inner and outer radii of the coaxial line, 
respectively.’ For electromagnetic analysis the geometry in 
Fig. 3(b) is equivalent to that in Fig. 3(a), [19]. Here a 
full-loop antenna is driven by a magnetic frill with the 
magnetic surface current density: 


M,(p) = -2V,/[pIn(a,/a;)]6, 


The electric field produced by this magnetic frill along the 
surface of the loop is the excitation for the antenna, £’. It is 
given by the following integral [20]: 


1 S 
-—|f M, x V’y dS’, 
4a S 


‘The radius of the outer conductor @, is chosen small enough that only 
the TEM mode propagates in the coaxial line. 


a;Sp=a,, 


a;<psa,. (8) 


Bi = 


(9) 


ZHOU AND SMITH: CIRCULAR HALF-LOOP ANTENNA 


where 


(10) 


and, referring to Fig. 4, r is the distance from the source 
point A (toroidal coordinates: p’, 6’, ¢’ = 0) on the frill to 
the observation point B (p = a,, 6, @) on the surface of the 
loop. The prime on V indicates differentiation with respect to 
the coordinates of the source point, and the surface of 
integration is the area of the frill. 

The field of the frill is most significant in the vicinity of the 
driving point, @ near ¢ = 0. Due to the thin-wire assump- 
tion, a; < b, the wire of the loop is nearly straight in this 
region, and the distance r can be approximated by 


r= y/4b? sin? (¢/2) + a? + p’? — 2a,p'cos(x) (11) 


with x = 6 — 6’. In addition, the component of the electric 
field tangential to the surface of the loop £%, is approximately 


(12) 


After inserting (8) into (9) and writing V’ in terms of the 
coordinates (p’, 6’, z’), E, is easily obtained, and (12) 


Ey, = E, cos ¢. 
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becomes 


— V, cos @ © 
BiG ea 
o(#) 2x In(a,/a;) | 


x 


tl 
| 
3 


ae V,cos 
2a In(a,/a;) 


f° 


x=-0 


[¥(o’ = @,) - ¥(0' = a;)] dx. (13) 


III. Fourrer SEr1ES SOLUTION 


The electric field of the excitation E4(o) can be expressed 
as a Fourier series: 


i Vo 
E,($) = s— 


es E toes b, cos(nd)}, — (14) 


with the coefficients given by 


b e 


b, = = 
. Vg g=-7 


E‘,(¢) cos (nd) do. (15) 


The electric field produced by the current (1) in the loop 
E%() (the scattered field) also can be expressed as a Fourier 
series. This has been done in the previous formulations 


[4)-[71: 


E3(o) = = [aol +2 e a,I,cos(n@)|, (16) 


where the coefficients a, are as in (4). When the boundary 
condition at the surface of the perfectly conducting loop is 
enforced: 


E4(o) + E3(¢) = 0, (17) 


and (14) and (16) are substituted, the Fourier series coeffi- 
cients for the current are obtained: 


-jV, [b 
ee (+). 
us a 


(18) 
n 

Notice that the Fourier series coefficients b, (15) for the 
delta-function excitation (2) are simply b, = 1, which with 
(18) leads to expression (3) for the Fourier series coefficients 
I, for the current. 

To obtain the current in the loop for the new method of 
excitation, the magnetic frill, the coefficients 5, first must be 
obtained. This requires evaluation of the integrals that result 
when (13) is substituted into (15) 


exp | —J6, V4? sin? (¢/2) + 4a? sin? (x /2) 


= sal [ 


—aY XS 8 


V 4B? sin? (¢/2) + 4a? sin? (x /2) 
exp | ~jB, 4b? sin’ (6/2) + a} + a? — 2a;a, cos (x) 


cos ¢ cos (nd) dx dd. (19) 


V4? sin? ($/2) + a? + a2 — 2a,a, cos (x) 


te a ag Sg vate B= 


qa Performance 
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The integration with respect to @ will be considered first; this 
can be split onto a number of integrals each of the form 


1 of 
ln e22% 
exp | —J8,b v/4sin? ($/2) + A?/b? + jpd| 


V4 sin? (¢/2) + A?/b? 


where A’/b? is a small quantity, such as 4(a;/b’) 
sin (x /2) < 1, and p is an integer. Wu evaluated similar 
integrals in obtaining the coefficients a, [4]-[7], and his 
result is 


K,( A) 


de (20) 


“2 0 [2,,,(x) + jJzp(x)] dx. (21) 


x= 


Wu includes the restriction ‘‘at least for p not too large’’ 
with the evaluation; however, (21) can be shown to hold for 
general p [18], [21]. 

Only the term that involves .%, in (21) contains the factor 
A. Thus when (21) is used with (19), the subtraction in the 
integrand cancels all terms except those involving %,: 


1 5 
co 2a In(a,/a;) i 


X= 


b 
x [AO Vina = (8.07 | 
Va? + a? — 2a;a, cos (x) 


x} : 


“V (ntl) - (6,6) | 
1] Va? + a, — 24,4, cos (x) 
b 


-V(n- 1) - 0% | dx. 


Now the integrals that remain can be evaluated using the 
following definite integral [22]: 


1 


27 


(22) 


KH,(c) £(d) = 


oO o 


af H,[ Ve? + d? — 2cdcos(t)] dt. (23) 


t=-3 
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The final result for the Fourier series coefficients for the 
electric field of the magnetic frill excitation is then 


ie wy (“lt Vint Ay = (6,6) | 


In(a, /a;) 


b 


aS Vn = 0 = (8,0) | 
- 4] V(n = 1) = (6,0) |}. 


Notice that these coefficients are no more complicated than 
the coefficients a, (4), (5); in fact, they are somewhat 
simpler since they do not involve the definite integrals of 0,,, 
and J,,,. 


+4] F vn 0? = (600) | 


(24) 


IV. NUMERICAL RESULTS 


With the Fourier series coefficients b, for the TEM coax- 
ial line excitation determined (24), the Fourier series coeffi- 
cients J, for the current in the half-loop (18) are readily 
obtained, and related quantities such as the input admittance 
(6) and radiated field [7] are easily calculated. 

In Fig. 2 the input admittance for the half-loop with TEM 
excitation is graphed as a function of the number of terms in 
the Fourier series. The parameters for the loop are as before, 
6,5 = 1.0, 2 = 10.0, and the dimensions of the coaxial line 
corresponds to a characteristic impedance Z, = 50 0, a,/a; 
= 2.30:+++ . The input conductance G of the loop is seen to 
be essentially the same as with the delta-function generator. 
However, the input susceptance B, unlike that for the delta- 
function generator, converges to a finite value as the number 
of terms in the Fourier series is increased. 

The convergence of the series for the input admittance can 
be examined by considering the limit as n— oo of the 
coefficients @, and 0, For a, in the limit n > oo: 


Ste 
Qa, = — P 
n Bb n 


with the following components of K,: 


285d 
| ‘ [Os.() + jJy,(x)] dx = 0, 
| 


7 +In [4(2? = p307)'"|} ea SF coer 


t 
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so, from (5): 


- tx, ).(2 1 aia 
n= — Ajo Nn Bi lied ——— 2 
T b 5” dam (=| 


thus, 


a, = — : 25 
278,48; (25) 


Similarly, for b,, in the limit 1 — oo 


ae waray 7) Fm) . x,( 5 n)| 


1 
In (a,/a;)(a,/b)n 


t 


or for n> b/(a, — a,), 
1 
er 
In (a,/a;)(a;/b)n 


The terms in the Fourier series for the input admittance (6), 
(18) are 


(26) 


—j2 { b, 
21,/V,= —|—}, n=l, (27) 
ot a, 
which for n > c becomes 
J46b 1 
aire tee (4). an 
53 In (a, /a;) n 


Now it is clear why the series for the input susceptance of 
the loop with TEM excitation converges: for large n it is 
proportional to 


Dae 


a series which is known to converge. It is also clear why the 
series for the input susceptance of the loop with the delta- 
function excitation does not converge. For the delta-function 
excitation, b, = 1 for all n, and (27) gives 


i4B,a; [1 
21,/%, = =| ). 


: (29) 


oO 


Thus for large n the series for the input susceptance is 
proportional to 


1 
eae 


a series which is known to diverge. 
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For practical calculations it is convenient to have an esti- 
mate of the error introduced when the Fourier series for the 
input admittance (6) is truncated at a finite number of terms 
N. Provided N > b/(a, — a;), from (28) this error will be 


e(N) = Y (infinite number of terms) — Y(N terms) 
—j2 Sad | b,, | J46,6 be 1 
a, f, In(a,/a;) n=N+1 ne 
(30) 


ous n=N+1 


To estimate the last summation, notice that 

oe. 1 ia 1 2N+1 1 
niNe1 W—-1 2N(N+1) N’ 
so 
J48,5 


(NYS ein(a, a) 


(31) 
The TEM field (7) is not the exact field in the aperture of 
the coaxial line; higher order modes are also present. For 
coaxial apertures of practical size, the higher order modes 
have a small effect on the input admittance of the antenna. 
The difference between the exact admittance Y and the 
admittance computed with the TEM excitation Yrpy,, 


5Y = Y— Yoem; (32) 


is smail and fairly independent of the type of wire antenna 
(linear, loop, etc.) excited by the aperture, provided the 
radius of curvature of the antenna wire near the aperture is 
large compared to the dimensions of the aperture. Morris 
determined 6Y for an infinitely long, linear antenna and 
presents his results in the form 


py el hot) (33) 
{ga 

where A is a constant depending on the ratio a;/a,, [23]. 
Values of 6 Y computed from his results are shown in Fig. 5. 
Notice that 6Y is small, a few tenths of a mS, for coaxial 
apertures of practical size (Z, ~ 50 9) and electrically thin 
wire antennas (8a; < 1). 

For the purpose of antenna design, it is convenient to have 
a graphical representation of the antenna input impedance 
Z=1/Y=R+/jX. In Fig. 6 the input impedance of the 
half-loop antenna is shown as a function of the electrical size 
of the loop 6,b = 27b/X,. Curves are given for different 
values of the thickness parameter Q = 2In(24b/a,). All 
results are for a driving coaxial line of characteristics 
impedance Z, = 50 Q and include the correction term 6Y 
(33); that is, 


Y = Yop + OY (34) 
is graphed. 


V. EXPERIMENTAL VERIFICATION 


An experimental model was constructed for the half-loop 
antenna, and measurements of the input admittance of the 
loop were made to verify the theory. The half-loop antenna 


ae oe oo ne Ma pert Es Bn eee 


Ba J. bau 
F. Manshadi 


1172 


Z,(0) 
138 72 41 31 6 


1.0 


0.8 


6,a,=0.150 
2) 
E 0.6 
mG 
> 
104 B,a,=0.075 
© 
0.2 
0.0 
0.0 0.2 0.4 0.6 0.8 1.0 
a/a, 
Fig. 5. Correction term 5Y versus the ratio of radii of inner and outer 


conductors of coaxial transmission line, a; /@,. 


was mounted on an aluminum image plane (155 cm x 120 
cm) and fed by a coaxial line from behind the plane. Absorb- 
ing material was placed around the front side of the image 
plane to reduce reflections from the support structure and 
edges of the plane. In addition, to further reduce errors due 
to reflections, the admittance was measured with the loop in 
eight different positions and these results averaged. 

The method of feeding the half-ioop is shown in Fig. 7. An 
APC-7 coaxial line feeds the antenna. A short spacer of 
length /, and radius a’, is placed in the end of the line to 
change the characteristic impedance of the feeding line from 
50 Q (APC-7) to Z,. The spacer and the loop wire have 
approximately the same radius. The spacer and loop wire 
were changed to obtain the three different values of character- 
istic impedance Z, and thickness parameter Q listed in Table 
I. 

The reflection coefficient for each loop was measured at a 
reference plane in the APC-7 coaxial line, Fig. 7. The 
reflection coefficient was transformed through the length of 
line / to the step discontinuity in the inner conductor. The 
admittance was then calculated and the capacitance of the 
step discontinuity removed. Values of this capacitance were 
calculated using the program of Jurkus [24] and are listed in 
Table I. Finally, to obtain an experimental admittance for 
comparison with theory, the corrected admittance was trans- 
formed through the length of the /, to the image plane. 

The measurements were made over the frequency range 
400 MHz < f<2 GHz using a Hewlett-Packard model 
8510A Network Analyzer with a Synthesized Sweeper. At 
each frequency point, 500 values of the reflection coefficient 
were measured and these averaged to reduce the effect of 
random noise. 

Results for the three loops listed in Table I are shown in 
Figs. 8(a)-8(c) where both G and B are graphed as a 
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(b) 
Input impedance of half-loop antenna as function of electrical size 
8,0. (a) Resistance (R). (b) Reactance (X). 


Fig. 6. 


function of the electrical size of the loop 6,b. The theory 
with the correction 6Y added (34) is shown as a solid line, 
while the experimental points are indicated by dots. The three 
cases correspond to feeding coaxial transmission lines of high 
(Z, ~ 86 ©), nominal (Z, ~ 46 Q), and low (Z, = 18 Q) 
characteristic impedance. In all cases the theoretical and 
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1, = 6.88mm 


Reference Plane 


Fig. 7. 
TABLE I 

PARAMETERS FOR THE HALF-Loop ANTENNAS USED IN THE EXPERIMENTS 

b 24; Q= 2a; a; Z, — Cap.of step 

Loop (cm) (mm) 21n(27b/a,;) (mm) a (Q) (pF) 

1 4.073 1.68 11.44 1.67 0.24 85.63 1.50 x 107? 
2 4.063 3.25 10.11 3.25 0.46 46.04 8.06 x 107? 
3 4.044 5.19 9.17 5.17 0.74 17.95 5.98 x 10-2 


experimental results are in excellent agreement. The slight 
difference that occurs in Fig. 8(c) for 8,b = 1.4 could be the 
result of the loop wire not meeting the thin-wire assumption 
of the theory. The thin-wire assumption requires B,a; = 
2xa;/, <1, the circumference of the loop wire must be 
small compared to the wavelength. At the highest frequency 
in Fig. 8(c) (8,5 = 1.7), the loop wire is over a tenth of a 
wavelength in circumference, 6,a; ~ 0.11. Some of the 
discrepancy for this case may also be due to the difficulty in 
fabricating such a fat loop so that it is circular and the 
difficulty in centering the fact loop/spacer in the coaxial line; 
the spacing between the inner and outer conductors was only 
0.9 mm. 

Figs. 9(a)-9(c) are circle diagrams for the input admit- 
tances of the three half-loops: B is plotted against G with 
G8, as a parameter. These graphs emphasize the behavior 
near the resonance of the loop, 8,b ~ 1.0, and again show 
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Image Plane 


Coaxial Line 
APC —7, Z, = 500 


3.04 mm 


2ay = 7.00mm 


Detail of method for feeding antenna in experiments. 


the excellent agreement between the theory and the experi- 
ment. 


VI. DiscussION AND CONCLUSION 


The conventional Fourier series analysis for the thin-wire 
circular transmitting loop, or its image equivalent the half- 
loop, uses a delta-function generator for excitation. The 
delta-function generator introduces two problems into the 
analysis: it does not correspond to any realizable method of 
feeding the antenna so an accurate comparison with measure- 
ment is not possible, and it produces a divergent series for 
the input susceptance B. A new model for the antenna was 
introduced that overcomes both of these difficulties: a half- 
loop driven through an image plane by a coaxial transmission 
line, with a TEM mode assumed in the aperture of the 
coaxial line. 

The Fourier series coefficients [, for the current in the 
loop now involve two factors: the coefficients a, (4) which 
appeared in the original analysis with the delta-function gen- 
erator, and the new coefficients b, (24); 


jeg (=). 
6m \a, 
Both a, and b, are easily calculated. They involve the 


modified Bessel functions 4, and .%, for which computer 
programs are readily available. In addition, @,, involves 
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Comparison of theoretical and measured input admittances (G, B) for half-loop antennas. (a) Loop 1, Q = 11.44; 


Zo = 85.63 Q. (b) Loop 2, 2 = 10.11; Zp = 46.04 Q. (c) Loop 3, 2 = 9.17; Zp = 17.95 Q. 
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Fig. 9. Circle diagrams comparing the theoretical and measured input admittances for half-loop antennas. (a) Loop 1, 2 = 11.44; 
Zo = 85.63 2. (b) Loop 2, 2 = 10.11; Zy = 46.04 Q. (c) Loop 3, @ = 9.17; Zp = 17.95 Q. 


definite integrals of the special functions Q,, and J,,, which 
may appear formidable; however, these integrals can be 
evaluated using the simple, rapidly convergent series pre- 
sented in the Appendix. 

The aforementioned (Introduction) extensions of the 
Fourier series analysis to other antenna structures can easily 
be modified to include the new method of excitation. For 
example, the full-loop antenna with concentric or eccentric 
spherical insulation and delta-function excitation has Fourier 
series coefficients for the current of the form [11], [13]: 


-JV, 


I, = ————_~., (35) The theory for the circular half-loop is now unique among 
Sm(a, + ¢,) the theories for thin-wire antennas; it provides simple, exper- 
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where c,, accounts for the insulation (this coefficient is gener- 
ally called b, in the previous analyses). The Fourier series 
coefficients for the half-loop antenna with hemispherical insu- 
lation and TEM excitation are simply 


a as (36) 


where b,, is given by (24). In (35) V, is the voltage of the 
delta-function generator; whereas, in (36) V, is the voltage at 
the end of the coaxial transmission line. 
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imentally verified, closed-form expressions (series) for the 
current and related quantities such as the input admittance 
and radiated field. 


APPENDIX 


The Lommel-Weber function has the series representation 
[25]: 


Dag(x) = = 
; ae Gre 8 ai 
[1? — (2n)’][3? - (2n)’] --- [(2k - 1) — (2n)7] 
(37) 


When this is substituted into the integral in (5) and integrated 
term by term, the following series is obtained: 


28,5 Lae 
[M0 de = - = 
W k=1 


, [-26,0)"|" = 
k[1? — (2n)’][3? - (2n)’] --- [(2k -— 1)" - Q22n)’] 
(38) 


Now this result can be expressed as 


28,6 tee 
[O° dal e)ae=-= OC, (39) 
x=0 W k=) 
where the recursive scheme 
— (28,b) (k — 1) 
k-1? (40) 


Gk - 17 = Gayl 


_ (28,6) 
2 eine (2n)’ 


(41) 


is used to compute the terms C,. 
A similar set of operations can be used to evaluate the 
integral in (5) of the Bessel function J,,: 
1 
T(2n + 1) 


x 


Jon(X) = (=)" 


(—x?/4)* 
kKi(2n+k)-+:(2n+1)] 


i¢> 


k=1 


(42) 


2(B,d)°"*' 1 


——— + 


Mone 
J,,( x) dx = 
x : 2n+1 k 


=0 T(2n+1) CK 


(43) 
where 


~(B,b)’(2n + 2k ~ 1) 


C, = = € 
« Qn+2k+ 1)k(Qn+k) *! (44) 
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1 


= ——., 45 
°  °€n4+1 (45) 


The series in (39) and (43) are convergent alternating 
series, and the errors made by neglecting all terms greater 
than the Nth are easily estimated: 


co 


bg 


k=) 


N 
Cy - es Cyl =|Cyail- 


Thus the absolute value of the relative error after summing 
N terms in either series is 


i) N 
Cee ee 
k=1 k=1 


N 
ae 
k=1 


|Cn-1| 


N 
Lc, 
k=1 


For loops of practical size, the series converge very rapidly; 
for example, with 6,b = 1.0 both series have a relative 
error less than 10~!° when nine terms are used. 
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